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Abstract 

We study the breaking of the string spanned between test charges in the three dimen- 
sional Abelian Higgs model with compact gauge field and fundamentally charged Higgs 
field at zero temperature. In agreement with current expectations we demonstrate that 
string breaking is associated with pairing of monopoles. However, the string breaking is 
not accompanied by an ordinary phase transition. 

1 Introduction 



> 

v~j 'The lattice Abelian Higgs model with compact gauge field (cAHM) in three dimensions is of 
•jjj .a broad interest both for high energy physics |2|] and condensed matter physics || |, H 
5^ ! where it was suggested to describe high-T c superconductors and strongly correlated electron 
systems. Nowadays, it has even entered the physics of cognitive networks [[]]. 

Due to compactness of the gauge field the model possesses Abelian monopoles which are 
instanton-like excitations in three space-time dimensions. The Abelian monopoles are able - if 
they are in the plasma state - to accomplish confinement of electrically charged particles. This 
is well known from cQED 3 where opposite charged particles are bound by a linear potential [[5] . 
The confinement is arranged by monopoles forming an opposite charged double sheet along 
the surface spanned by the trajectories of the external test charges. This surface is usually 
considered as the world surface of a string. Due to screening, the free energy increases only 
proportional to the area of the surface such that an area law for the Wilson loop emerges. 

However, if dynamical matter fields in the same representation as the external test charges 
are added to the confining theory, linear confinement may be lost. This should be so, irre- 
spective whether the dynamical matter field is fermionic (the quarks in QCD) or bosonic (the 
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Higgs particle in our case). The string breaking phenomenon has been extensively studied in 
non-Abelian gauge theories with matter fields [§] or with test charges in the adjoint represen- 
tation ||. Here we want to investigate string breaking in cAHM 3 with a q = 1 charged Higgs 
field, a model whose permanently confining counterpart, cQED 3 , is well understood. The gen- 
eral, intuitive picture says in the present case that the string breaks because of Higgs particle 
pairs popping up out of the vacuum at a definite inter-particle separation between the external, 
infinitely heavy test charges. Thus, the physical state corresponding to a broken string would 
consist of two heavy-light mesonic states plus some number of light-light Higgs pairs. 

In order to destroy the linearly rising potential in cAHM 3 , the coupling between the Higgs 
field and the gauge field must be sufficiently strong. One might be tempted to associate the 
string breaking with a phase transition between confinement and Higgs phases. Indeed, Ref. || 
proposes to associate the string breaking with a Berezinsky-Kosterlitz-Thouless type transition. 
In this paper we demonstrate that and how the expected string breaking happens in a part of 
the phase diagram where a first or second order phase transition can definitely be excluded. 
Abelian monopoles play the crucial role in the dual superconductivity scenario [XD| of con- 



finement in QCD. There, the monopole degrees of freedom need to be defined with the help of 
Abelian projections ]TT[ (see, e.g. reviews The condensed magnetic currents were shown 

to make a dominant contribution to the string tension between quarks, in pure SU(2) gauge 
theory [JIB as well as in SU(3) gluodynamics and also in full finite-temperature QCD with 
Nf = 2 flavors of dynamical quarks fl4[]. Moreover, in full QCD with dynamical quarks the 



contribution of Abelian monopoles to the heavy-quark potential QCD shows the property of 



string breaking |T5|. The breaking of the adjoint string in pure gluodynamics as well as the 
breaking of the fundamental string in full QCD can both be described within the Abelian pro- 
jection formalism [[L£|. The back-reaction of the dynamical fermions on the gauge field should 
modify the dynamics of monopoles in such a way that this dynamics incorporates the above 
qualitative picture ]15| . 



Therefore, guided by the analogy to QCD, we focus our interest in the present paper on 
the monopole degrees of freedom in compact AHM in three dimensions under the influence of 
a scalar matter field. We would like to elucidate the changing role of monopoles under the 
particular aspect of string breaking. As in QCD, the string tension in this model is exclusively 
due to monopoles. Therefore one can expect that monopoles also encode the back-reaction 
of the matter field causing the string breaking phenomenon. Here we want to demonstrate 
that (i) the monopole part of the potential indeed incorporates the effect of string breaking 
and (ii) that it is monopole pairing which is the reason for the breakdown of the monopole 
confinement mechanism. We are aware of the incompleteness of the analogy to QCD and the 
relative simplicity of monopole dynamics in 3 instead of 4 dimensions. 

It seems that there is only one possibility to explain string breaking in three space-time 
dimensions. We assume that, in the presence of matter fields, monopoles are increasingly bound 
into neutral pairs (magnetic dipoles). The size of a typical pair should be of the order of the 
string breaking distance R^ r . Indeed, if the distance R between the test charges is much larger 
than i?br then the test charges do not recognize individual monopoles inside the dipoles (in other 
words, the fields of the monopoles from the same magnetic dipole effectively screen each other) 
and the vacuum is basically composed of neutral particles. Therefore, at large inter-particle 
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separations there should be no string tension. However, if R <C -Rbr then the test charges do 
recognize individual monopoles even if they are bound in dipoles, and the monopole fields may 
induce a piecewise linearly rising potential. These simple considerations can be made more 
rigorous by analytical calculations |L7| for a gas of infinitely small-sized dipoles. 

Recently, it was found that the matter fields in the Abelian Higgs model lead to a loga- 
rithmic attraction between monopoles and anti-monopoles 0] which results in the formation 
of monopole-anti-monopole bound states and string breaking. The formation of dipoles can 
also be explained as due to the existence of Abrikosov-Nielsen-Olesen vortices | i~8f , the string 



tension of which gets increased as we move in the parameter space deeper into the Higgs re- 
gion []. Massless quarks also force the Abelian monopoles to form bound states |1| . Note that 
the origin of monopole binding in the zero temperature case of CAHM3 is physically different 
from the monopole binding observed at the finite temperature phase transition in compact 
QED pifl . It is different as well from the Z2 vortex mechanism in the Georgi-Glashow 
model J22[. 



In this paper we numerically establish a relation between string breaking on one hand and the 
occurrence of monopole-antimonopole bound states on the other by studying some properties 
of the monopole ensembles provided by the compact Abelian Higgs model. In Section [2] we 
recall the definition of the model and discuss its missing ordinary phase transition. In Section |3] 
flattening of the potential is described. Here we also introduce the i] angle as a parameter which 
defines the "effectiveness" of string breaking. Section f| is devoted to an investigation of the 
cluster structure of the monopole ensembles. Our conclusions are presented in the last Section. 



2 The Model and Its Crossover 

We consider the 3D Abelian gauge model with a compact gauge field 9 Xtfl and a Higgs field 
§ x with unit electric charge. The coupling between the gauge and the Higgs fields is S X)li oc 
3?e (§ x e iex ^§ x+ ( l ). To simplify calculations we consider the London limit of the model, which 
corresponds to an infinitely deep potential on the Higgs field. In this limit the radial part 
of the Higgs field, \<& x \, is frozen and the only dynamical variable is the phase ip x of this 
field, <& x = \$ x \ e lipx . Thus the Higgs-gauge coupling reduces to the simple interaction S Xtfl oc 
cos(ip x+ f l — ip x + 9 X4 j). However, the model can be simplified even further by fixing the unitary 
gauge, ip x = leading to S x ^ cx cos 6^. Thus we consider the model with the action 

3[9] = -(3j2 cose P- K 12 C08d i> W 
p 1 

where (3 is the gauge (Wilson) coupling, k is the hopping parameter and dp is the plaquette 
angle. We study the model at zero temperatures on lattices of size L 3 , with L = 12, 16, 24, 32. 

The phase structure of the model on the boundaries of the phase diagram in the /3— k 
plane can be established using the following simple arguments. At zero value of the hopping 
parameter k the model ([l]) reduces to the pure compact Abelian gauge theory which is known to 

^^Note that the division of the parameter space of the model into Higgs and confinement regions is only loose 
since these regions - as we discuss below - are analytically connected. 
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be confining at any coupling f3 due to the presence of the monopole plasma |7J]. This argument 
extends to the \ow-k region of the phase diagram. Therefore we call this the "confinement 
region" . At large values of k (also called the " Higgs region" ) the monopoles should disappear 
because the gauge field in this limit is increasingly restricted to the trivial vacuum state: 



8 x ,n — 0. 



At large /3 the model reduces to the three dimensional XY model which is known to have 
a second order phase transition at k* y « 0.453 f23 |. Indeed, in this limit we get the condition 



d9i = 6p = which forces the gauge field to be a gauge transformation of the vacuum, 9 Xtfl = 
—<Px+fi + <Px + 27r/ XjAt G (—it, 7r], l X)lt G X, (p x G (— 7r, tt]. The scalar fields are the spin fields in 
that model. 

Despite the phase structure on the boundary of the coupling plane is well established, the 
structure of its interior is still under debate. Indeed, in Ref. [[| arguments were given that 
the interior is trivial (i.e., there is no ordinary phase transition for finite values of /3 and k) 



while the XF-phase transition takes place in an isolated point at j3 — oo. In Ref. it has 
been suggested that the phase diagram of cAHM 3 resembles the vapor-liquid diagram with 
a critical end-point. Finally, in Ref. [0] it was argued that the phase diagram contains a 
"pocket" in which a Coulomb phase could be realized. Arguments given in Ref. || do not allow 
to distinguish between these three possibilities. 



In a numerical study on rather small lattices G9] no hint for an ordinary phase transition 



at finite coupling constant f3 has been found. However, for simulations allowing fluctuating 
Higgs lengths, sufficiently away from the London limit, the phase diagram has been seen to 



become nontrivial [pq| . Recently, the phase structure of the CAHM3 has been studied by the 
authors of Ref. [27] in connection with the nature of the transition in the type-I and the type-II 
region. The alleged second order transition in the type-II region away from the London limit 
still remained inconclusive. 

Here we are not going to study the whole phase diagram of cAHM 3 although this question 
would be still interesting. As we describe below, we observed that at moderately small (3 
the Higgs and confinement regions are connected analytically by a crossover as predicted in 
Ref. 0. We concentrate on the changing role of monopoles under the aspect of the string 
breaking phenomenon accompanying the crossover at relatively small (3 with increasing hopping 
parameter k. For the simulations we use a Monte Carlo algorithm similar to the one described 
in Ref. JHJ and have considered 5 ■ 10 3 to 5 ■ 10 4 independent configurations per data point, 
depending on the lattice size and the set of coupling constants. We vary the value of the hopping 
parameter k at a fixed value of gauge coupling constant f3 = 2.0. To locate a (pseudo-) critical 
point we use the susceptibility of the hopping term, 

X=(S 2 H [9]}-(S H [9}}\ SW[0] = -£cos0,, (2) 

1 

which is shown in Figure |l](a) for L = 12, 16, 24, 32. The height of the peak is practically 
independent on the lattice size. We have observed a very similar volume independence also of 
the susceptibility of the gauge term, Sg[0] = —J2p cos Qp- Thus, in agreement with Ref. |25j we 



2 Note that all figures in this paper are shown for (3 = 2.0. 
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conclude that there is no ordinary phase transition between the Higgs and confinement regions 
of the parameter space of the model. 

The crossover point k c (L) is located fitting the susceptibility (§) in the vicinity of the peak 
by the following function: 



Ci 



[ci + (k- n c y\ 



2-i Q ' 



(3) 



where Ci^, k c and the power a are fitting parameters. In Figure |T](a) we show the fit of the 
susceptibility data for the 32 2 lattice. The fit parameters practically do not depend on the 
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Figure 1: (a) The susceptibility of the hopping term (0) as a function of k; (b) the crossover 
point function of the inverse lattice size, L . 

lattice size. We depict the critical value of the hopping parameter k c vs. the inverse lattice size 
L^ 1 in Figure |l|(b). The value of the power a is very close to 1/4. In the next Sections we will 
work with the lattice 32 2 which passes the crossover at k c = 0.526(1) along a line with fixed 
P = 2.0. 



3 The Flattening of the Potential 

String breaking manifests itself in the flattening of the potential between test particles with 
(opposite) electric charges q — ±1. In principle, we can separate the contributions to the 
potential from monopoles and from the rest ("photon contribution"). Monopoles are responsible 
for the string tension. Therefore one can expect that the monopole contribution alone will signal 
the onset of string breaking when the monopole dynamics starts changing. It would be much 
more demanding to extract the string part from the full potential and to study its change over 
the parameter space of the model. The full potential contains also the perturbative photon 
contribution which - being logarithmically large at small distances - shadows the eventually 
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linearly rising part. Any statement about the string part would require a careful fit of full 
potential. On the more technical side, the monopole contribution alone, calculated separately 
according to the configurations generated in the simulation of the AHM, has a much better 
signal/noise ratio compared to the full potential. All this justifies to proceed directly to the 
evaluation of the monopole contributions to the external-charge potential. 

To this end we have divided the gauge field 9 t into a regular (photon) part and a singular 
(monopole) part 



e = e phot + e mon , e mon = 2n^ l 5p[i] . (4) 

The 0-form *j G Zj is nonvanishing on the sites dual to the lattice cubes c which are occupied 



by monopoles ]29 



3c = ±z £ (-ll'^U. (5) 

where the factor (— 1) P takes the plaquette orientations relative to the boundary of the cube 
into account. In Eq. @) the 2-form pp[j] = [Bp] (the notation [• ■ •] means taking the integer 
part) corresponds to the Dirac strings living on the links of the dual lattice, which are either 
closed or connecting monopoles with anti-monopoles, = *j. While *j is gauge invariant, 

the 2-form pp[j] is not. For the Monte Carlo configurations provided by the simulation of ([!]) 
we have located the Dirac strings, p[j] 7^ 0, and constructed the monopole part mon of the 
gauge field according to the last equation in (^). The operator A3 1 in Eq. (|j) is the inverse 
lattice Laplacian defined for a three-dimensional lattice L 3 : 

A7 1 (x;L) = — -Y —, , (6) 

where p, = 2-nkijLi for ki = 0, . . . , Li — 1, with % — 1, . . . , d and Li = L. 

We define the potential between test particles with the help of the following correlator of 
two Polyakov loops: 

{P(6)P\R)} = e~ LViR) , (7) 

located at two-dimensional points and R. The potential V depends on R = \R\. The 
use the Polyakov loop has clear advantages compared to the Wilson loops. The construction 
of the Polyakov loops is not only possible for finite-temperature but also for finite-volume 
cases. L = Li is the common length of the zero-temperature box in all three directions. Due 
to the absence of space-like links joining the Polyakov loops the correlator ([FD defines the 
static component of the potential. Note that the monopole contribution to the Polyakov loop 
correlator (0) does not depend on the precise form of the Dirac string *p[j}. Therefore this 
contribution is gauge-invariant. 

We discuss the results for the potential using the following fitting function: 



e 



^ = Co [siu 2 r, + cos 2 r, C ° Ma l\ L 'l~ R)) ] -exp^A^-A^O)]}, (8) 

cosh(cr///2) J 1 - " - ' 
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where Co, T], a and 7 are fitting parameters and A2 is the inverse lattice Laplacian in two 
dimensions. 

The meaning of the expression (H) is quite simple. In the absence of string breaking and 
in an infinite two-dimensional volume the leading contribution to the function in the right 
hand side of Eq. @ should be just const • e~ aLR where a is the effective string tension. Here 
"effective" means that this term gives rise to a linear part in the potential at short distances. 

The string breaking manifests itself in the appearance of an additional constant term, 
consti + const?, ■ e~ aLR . Next, the finiteness of the two-dimensional volume reduces the ex- 
ponential to the cosh-function which takes care of the symmetry R — > L — R. Finally, we 
introduced a Coulomb term in order to take into account sub-leading corrections. 

The dimensionless parameter 77 e [0, tt/2] - which we call a "breaking angle" - has a sense 
only as long as a 7^ 0. It can be considered as a kind of "order parameter" for string breaking: 
if r\ — 0, no string breaking occurs, and if rj = 7r/2, the potential does not contain a linear piece 
at all. An intermediate value of the breaking angle implies the existence of the finite distance 
i? s b at which the string between the test particles breaks. Note that we have introduced a 
normalizing cosh-factor in the second term in the brackets in order to keep the V Rt (R = 0) 
value independent on 77. This definition is a matter of conventions. 

To justify the presence of the Coulomb-like term in the fitting function let us con- 
sider three dimensional compact QED. It is well known that in the Villain representation the 
Polyakov loop correlator factorizes into the photon and monopole contribution. The monopole 
contribution can be evaluated exactly and it contains a massless pole, A^" 1 (i?), correspond- 
ing to the Coulomb potential between test particles. The total correlator should not contain 
the massless pole due to the massiveness of the photon. Therefore the monopole contribution 
to the correlator must contain - in addition to the linear term - the difference between the 
Yukawa and Coulomb potentials, A? 1 (R;m) — A? 1 (R;m = 0) corresponding to the exchange 
by "real" (massive) and "bare" (massless) photons. Here A% (R;m) is the propagator of a 
particle with the mass m. The mentioned above sub-leading term is small at distances smaller 
than the inverse photon mass. However, this term gives a significant (logarithmically growing) 
contribution at larger separations between test particles. Thus the largest deviation from the 
linear behaviour of the monopole contribution to the potential is expected to come from large 
distances due to exchange of a massless (bare) photon. 

Similar arguments should apply to the case of the compact AHM. The bare photon here, 
however, is not massless due to the spontaneous breaking of the U(l) symmetry. Therefore 
the fitting function (^j) should be modified: the Coulomb potential should be replaced by the 
Yukawa one. We have found that such fits do not work well because the corresponding Yukawa 
mass turns out to be consistent with zero within huge error bars. On the other hand, the mass 
of the bare photon should be small at the QED side of the crossover where the form of the fit (^) 
is obviously justified. We have found numerically that this fitting function works well also at 
the Higgs side of the crossover. Therefore in Eq. (§) we restrict ourselves to the Coulomb term 
only. 

The fits of the numerical data for the potential V(R) due to monopoles by the expression 
(§) are shown in Figure ^|(a) for five values of the hopping parameter from k = 0.52 (below 
string breaking) to k = 0.60 (far from the transition on the Higgs side) including k = 0.53 m k c 
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Figure 2: (a) The potential for k = 0.52, k = 0.53, k = 0.54, k = 0.55 and k = 0.60 extracted 
from the monopole contribution to the Polyakov loop correlator by Eq. ([?]). The fits by the 
function (§) are shown by solid lines, (b) The string tension vs. k. In this and all subsequent 
figures the string breaking transition at k c for f3 = 2.0 is marked by a vertical line. 

(in the vicinity of the transition). In the fits of the potential the point R = was excluded. 
One can clearly recognize a linear part in the potential near the transition point. As k increases 
(this corresponds to moving deeper into the Higgs region) the linear part gradually disappears. 
This can also be seen from the properties of the string tension a shown in Figure |2](b). The 
string tension itself, which on the confinement side amounts roughly to 50 % of the QED3 string 
tension (corresponding to k = 0), drops to a smaller value over a very narrow k region. The 
described behaviour of the potential is consistent with the expected disappearance of isolated 
monopoles on the Higgs side of the string breaking transition. The residual string tension, 
which is accompanied by a short string breaking length i? s b, can be accounted for by the 
monopole-antimonopole dipoles of finite size. With 77 — 7r/2 the fit error of a increases. 

The breaking angle r] is shown in Figure |3|(a) as a function of k. It clearly shows an "order- 
parameter-like" behaviour: it is close to zero for k < k c and it is finite at k > k c . Small 
values of rj imply that the string breaking distance is still large. At k ~ 1 the value of 77 ~ 7r/2 
indicates that the area-law term in the Polyakov loop correlator (|S]) has become irrelevant. 

The parameter 7, shown in Figure 0(b), seems to vanish on the Higgs side of the string 
breaking transition. This may indicate that in the Higgs region the "bare" photon mass becomes 
significant and that the corrections to the linear potential gets concentrated at small distances 0. 
Thus, long distance corrections should be zero, i.e. 7 ~ 0. 

3 We remind the reader that the smallest distance, R = 0, is excluded from the fit. 



8 



jt/2 
X] 

jt/4 



-0.05 



0.4 K 0.6 
c 



0.8 K 1 



(a) 




-0.1 



Figure 3: (a) The breaking angle i] and (b) the parameter 7 appearing in the fitting function (||) 
vs. the hopping parameter k. 

4 The Cluster Structure of the Monopole Ensemble 

In this section we turn to the monopole clustering aspect of the Monte-Carlo configurations 
which have been used in the last section to work out the monopole part of the external-charge 
potential. We closely follow Ref. |ZT| where the cluster analysis of the monopole configurations 
in the case of compact QED3 at non-zero temperature was performed. 

The simplest quantity describing the behaviour of the monopoles is the monopole density, 
P = J2 c \jc\/L 3 , where j c is the integer valued monopole charge inside the cube c defined in 
Eq. flU). The density of the total number of monopoles is a decreasing function of the hopping 
parameter k as it is shown in Figure f|(a) by diamonds. The density sharply drops down at 
k c , which has been recognized as the string breaking transition point, but the density does 
not vanish on the Higgs side of the crossover. The binding of monopoles into dipoles should 
show up as an increase of the number of monopoles enclosed in neutral clusters. We call a 
monopole cluster neutral if the charges of the corresponding constituent monopoles sum up to 
zero. Clusters are connected groups of monopoles and anti-monopoles where each object is 
separated from at least one neighbor belonging to the same cluster by a distance less or equal 
than some R ma , x . The smallest clusters are isolated (anti-)monopoles. In our analysis we have 
used -R^ax = 3 a 2 which means that monopoles are considered as neighbors if their cubes share 
at least one single corner. 

We show also in Figure ||(a), symbolized by triangles, the density of monopoles in neutral 
clusters which almost covers the total density on the Higgs side of the string breaking transition. 
If we take into account that also bigger dipoles - which cannot be identified by our procedure 
- may be formed, this clearly signals the binding transition. 

In an alternative, perhaps more clear way this is illustrated by the fraction of monopoles 
belonging to neutral clusters, N = p n0 utrai/ptotai, which is shown in Figure f|(b). Being constant 
on the confinement side of the string breaking transition, this quantity starts suddenly to rise 
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Figure 4: (a) Total density of monopoles (diamonds) and density of the monopoles in neutral 
clusters (triangles), and (b) the fraction of neutral clusters among all clusters, both as functions 
of K. 



at the transition. This indicates that at the transition point (crossover) the binding process 
rapidly takes place. At large k the fraction is very close to unity. Then all monopoles are 
bound. 
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Figure 5: (a) The average number of monopoles and ant i- monopoles per cluster as function of 
k and (b) the (normalized) cluster size distribution D(s) for a few values of k. 



Finally, in Figure |5](a) we present the average number of (anti-)monopoles per cluster and 
in Figure ||(b) the (normalized) cluster size distribution D(s) where s is the number of (anti-) 
monopoles in the cluster, for a few values of the hopping parameter k. On the confinement side 
of the string breaking transition (re pa 0.5) the vacuum consists to m 70 % of isolated monopoles. 



10 



At the crossover (the string breaking transition) at k ~ 0.53 the number of isolated monopoles 
decreases, and on the Higgs side (k > 0.53) the vacuum is dominated by the dipole gas. 



5 Conclusions 

We have numerically observed that in the London limit of the three-dimensional Abelian Higgs 
model string breaking occurs and is accompanied by monopole recombination into dipoles, in 
agreement with arguments given in Ref. [|J. 

Our study shows that the monopole binding is not necessarily accompanied by an ordinary 
phase transition of first or second order. There is a proposition [|j], however, that the string 
breaking may be associated with a Berezinsky-Kosterlitz-Thouless type transition due to the 
appearance of an anomalous dimension of the gauge field induced by the fluctuations of the 
matter fields. This possibility is not ruled out by our results. In the London limit (studied in 
this article) the fluctuations of the radial components of the matter field are suppressed, while 
far away from the London limit the fluctuations become significant such that an ordinary phase 
transition may exist [F2^, CT . 
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